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Stability  and  normal  zone  propagation  speed  in 
YBCO  coated  conductors  with  increased 
interfacial  resistance. 


George  A.  Levin,  Paul  N.  Barnes,  Jose  P.  Rodriguez,  Jake  A.  Connors,  and  John  S.  Bulmer 


Abstract —  We  will  discuss  how  stability  and  speed  of  normal 
zone  propagation  in  YBCO-coated  conductors  is  affected  by 
interfacial  resistance  between  the  superconducting  film  and  the 
stabilizer.  Our  numerical  simulation  has  shown  that  the 
increased  interfacial  resistance  substantially  increases  speed  of 
normal  zone  propagation  and  decreases  the  stability  margins. 
Optimization  of  the  value  of  the  resistance  may  lead  to  a  better 
compromise  between  stability  and  quench  protection 
requirements  than  what  is  found  in  currently  manufactured 
coated  conductors. 

Index  Terms — Coated  Conductors,  Normal  Zone,  Stability. 

I.  INTRODUCTION 

uperconducting  wires  in  large  scale  applications  have  to 
satisfy  two  important  criteria  -  they  have  to  be  stable 
during  operation,  but  also  allow  integration  of  an 
effective  quench  protection  scheme  in  the  superconducting 
devices  [1].  These  two  requirements  compete  with  each  other. 
Conventional  low  temperature  superconductors  operate  at 
liquid  helium  temperatures,  where  their  small  heat  capacity 
makes  the  current-carrying  wires  very  unstable  with  respect  to 
even  a  small  amount  of  heat  released  in  them.  For  that  very 
reason,  however,  the  speed  of  normal  zone  propagation  (NZP) 
is  high  (several  meters  per  second)  which  is  beneficial  for  the 
rapid  detection  of  the  normal  zone  (NZ)  and  mitigation  of  the 
consequences  of  the  quench.  The  situation  is  reversed  in 
YBa2Cu307.x  (YBCO)  coated  conductors  [2].  Their  main 
advantage  over  conventional  superconductors  is  a  high 
operating  temperature  (60-77  K)  and  the  resultant  high 
stability  margins.  The  flip  side  is  a  low  NZP  speed  which 
complicates  quench  protection  [3-5]. 

We  present  the  results  of  a  numerical  analysis  of  a  model 
of  NZP  specialized  to  the  architecture  of  the  state  of  the  art 
coated  conductors.  The  main  conclusion  of  this  analysis  is  that 
the  NZP  speed  can  be  increased  up  to  an  order  of  magnitude 
by  increasing  the  interfacial  resistance  between  copper 
stabilizer  and  the  superconducting  film.  The  concomitant 
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reduction  of  the  stability  margin  is  substantial,  but  the  wire 
retains  an  ability  to  recover  on  its  own  from  a  finite  amount  of 
locally  injected  heat. 

After  some  of  our  preliminary  results  were  presented 
elsewhere  [6,7],  we  were  advised  by  colleagues  that  a  large 
body  of  work  devoted  to  the  effect  of  interfacial  resistance  on 
NZP  in  conventional  superconductors  had  been  done  much 
earlier,  with  some  of  the  publications  predating  the  discovery 
of  high  temperature  superconductivity  [8-11].  This  reminds  us 
again  of  an  ancient  adage  that  -e-ach  new  thing  is  a  well- 
forgotten  old  one”. 

II.  MODEL 

YBCO  coated  conductors  are  thin  metal  tapes  in  which  an 
about  1  pm  thick  YBCO  film  is  sandwiched  between  a  metal 
substrate  and  the  copper  stabilizer  [2,12].  The  total  thickness 
of  the  tape  is  about  evenly  divided  between  stabilizer  and 
substrate  and  close  to  100  pm.  The  standard  width  of  the 
currently  manufactured  tape  is  about  4  mm.  There  is  a  very 
thin  layer  of  material  between  the  YBCO  and  copper  that 
accounts  for  the  resistance  to  the  current  exchange  [13]. 
Common  wisdom  is  that  the  resistance  of  this  interface  has  to 
be  as  low  as  possible  to  ensure  effective  current  transfer.  If 
necessary,  however,  it  can  be  increased  by  a  chemical 
treatment  or  by  changing  the  architecture  of  the  wire.  Here  we 
treat  the  interface  as  an  infinitesimally  thin  boundary  between 
YBCO  and  copper  with  a  finite  resistance  R  [Q  cm^]. 

Hereafter,  we  consider  only  the  direct  transport  current 
condition.  Three-dimensional  (3D)  equations  of  heat 
conduction  in  a  thin  tape-like  composite  wire.  Fig.  1,  can  be 
reduced  to  the  2D  (planar)  or  ID  linear  model  if  the  variation 
of  temperature  along  the  thickness  of  the  wire  is  small  in 
comparison  with  its  variation  along  the  plane  of  the  wire.  One 
can  show  that  this  criterion  is  satisfied  if 
djr  «K.{Tc—Tq)  ,  where  yV  is  the  thermal  flux  lifted 

from  the  surface  of  conductor,  and  dt  and  Kf,  respectively,  are 
the  thickness  and  thermal  conductivity  of  either  the  stabilizer 
or  substrate  material;  Tc  and  Tq  are  the  critical  temperature  and 
the  operating  temperature  respectively  {Tc  -Tq  ^20  K).  For 
copper  stabilizer  Ki  »  4-5  W/cm  K  and  di  ^  40  pm.  For  a 
substrate  metal  like  Hastelloy  K2  »  7x10'^  W/cm  K  and 
typically  ^2  ^  50-100  pm  [13].  As  long  as  the  heat  flux  from 
the  surface  is  of  the  order  of  1-10  W/cm^  or  less,  the  planar 
approximation  is  justified.  The  150-300  nm  thick  buffer 
between  the  YBCO  and  substrate  and  even  thinner  YBCO-Cu 
interface  also  can  be  considered  thermally  transparent. 
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The  2D  -  {x,y}  -  heat  conduction  equation  can  be  obtained 
by  integrating  the  3D  equation  along  the  thickness  of  the 
conductor  (Fig.  1): 

c  km  =  Q-2K,{T-T,)  (1) 

ot 

Here  C=C]d]+C2d2  is  the  combined  specific  heat, 
K=Kidi+K2d2  is  the  effective  thermal  conductivity,  and 
Q=fq(z)dz  is  the  areal  density  of  the  internal  heat  sources.  Kg 
is  the  heat  transfer  coefficient  across  the  insulation  on  the 
surface  and  Tg  is  the  ambient  temperature  which,  in  the 
absence  of  losses,  is  the  operating  temperature.  We  also 
neglect  here  for  simplicity  the  temperature  dependence  of  the 
material  coefficients. 

To  obtain  a  closed  system  of  equations  for  a  simplified 
analysis  we  shall  neglect  the  dissipation  {Q  =  0)  in  the  sections 
of  the  wire  that  are  subcritical,  J  <  J/T).  This  is  an 
approximation  because  a  certain  amount  of  current  flows  in 
the  stabilizer  even  in  the  subcritical  sections  [12].  The  limits 
of  this  approximation  are  discussed  below.  We  shall  employ 
the  Bean  model  approximation  to  describe  the  internal  heat 
sources  in  the  sections  of  the  wire  that  are  still 
superconducting,  but  where  the  transport  current  density 
exceeds  the  local  critical  current  density,  J  >  MT)  [9], 
Hereafter  the  symbols  J  denote  the  sheet  current  density 
[A/cm  width].  There  are  three  separate  contributions 
originating  from  the  stabilizer,  the  interface  and  the 
superconductor,  respectively: 

q=—P+  Siz)  +  J, \E^  |^(z) .  (2) 

A  R 

Here  is  the  electric  field  in  the  stabilizer  and  pi  is  its 
resistivity.  Vi  is  the  local  potential  of  the  stabilizer  and  Vs  is 
the  potential  of  the  superconductor.  The  5-functions  account 
for  the  fact  that  two  of  the  heat  sources  are  concentrated  in  the 
volume  much  thinner  than  either  the  stabilizer  or  substrate. 
The  heat  dissipation  in  the  substrate  is  negligible  due  to  its 
large  resistance  [12]  and  the  substrate  contributes  only  to  the 
thermal  mass  of  the  conductor  (especially  in  light  of  the  thin 
buffer  layers). 

Within  the  Bean  model  approximation  the  current  through 
the  superconducting  film  cannot  exceed  the  critical  current,  so 

that  Js=5Jc(T)  ,  where  n  is  a  unit  vector  in  the  direction  of 

current.  The  excess  of  transport  current  must  flow  into  the 
stabilizer.  The  condition  of  charge 

conservation  V  •  =  0 ,  where  -  VJ/R  is  the 

density  of  current  across  the  interface,  takes  form 

«-VJ,=(l^-VJ/R  (3) 

In  the  one-dimensional  case  that  will  be  considered  hereafter, 
Eq.  (3)  and  its  spatial  derivative  are  given  by 

//=(F,-VJ/R;  RJf  (4) 

Here// /ax  and  j'^=d^JJdx^‘  The  direction  of  the 

current  flow  is  taken  to  be  positive.  Fig.  1.  Therefore,  Eg  must 
be  either  positive  or  zero. 

The  power  sources  in  the  right-hand  side  of  Eq.  (1)  have  a 
different  form  in  the  three  temperature  regions:  T<Ti, 


Tj<T<Tc  ,  and  T  >  T^,  respectively.  The  current  sharing 
temperature  T/  is  defined  by  the  condition  Jc(T i)=J.  For  T<Tj 
(sub-critical  regime)  we  take  Q=0.  We  want  to  emphasize  that 
this  is  an  approximation  we  have  mentioned  earlier.  In  fact 
when  resistance  of  the  interface  is  large,  the  current  flows  in 
the  stabilizer  even  in  the  sufficiently  -eold”  sections  where 
Jc(T)  >  J.  The  critical  temperature  7/  is  defined  by  the 
condition  Jc(Tc)=0.  At  temperatures  above  7^  the  transport 
current  flows  only  in  the  stabilizer  and  Q=pf/di.  All  three 
heat  sources  in  Eq.  (2)  contribute  in  the  intermediate  range  of 
temperature,  Ti<T<Tc. 


Q  = 


dy 


+  R{Jlf+J, 


d^ 


H{E^) 


.(5) 

A  unit  step  function  h{E  )  ensures  that  the  electric  field  in  the 

superconductor  is  either  positive  or  zero.  The  necessity  to 
introduce  the  step  function  -by  hand”  stems  from  the 
approximation  mentioned  above,  when  we  neglected  the 
precursor  electric  field  propagating  ahead  of  the  normal  zone 
[12]  in  the  sections  of  wire  where  7<7/.  The  step  function 
ensures  that  the  numerical  solution  does  not  zander  off’  into 
unphysical  territory  where  the  heat  source  may  become 
negative. 

Rewriting  Eq.  (1)  in  dimensionless  form  will  allow  us  to 
combine  numerous  material  and  experimental  parameters  into 
a  smaller  number  of  dimensionless  quantities  that  determine 
different  regimes  of  the  NZP.  A  dimensionless  temperature  0 
is  defined  as 


0  J^{0)  =  J{\-9).  (6) 

E-T, 

Here  we  also  assumed  for  simplicity  a  linear  dependence  of 
the  critical  current  on  temperature,  well  justified  for  coated 
conductors  [13].  The  distances  shall  be  measured  in  units  of 
the  thermal  diffusion  length  =  (7)^  /  ?  where 

—  K I C ,  and  time  in  units  of  y“\  where  the  increment 
Y  —  !  d^Chl  determines  the  characteristic  time 

required  to  warm  an  element  of  the  conductor  by  A7  -T^-T^. 


Equation  (1)  takes  the  following  piece-wise  form: 


de 


d^9 


dr  df  “ 


d>l. 


(7) 


--^  =  0^+  r{9' )"+(!-  0)(e  -  r0"  )H{9  -  r0"  )-k(0-  0, ); 
dz 

0<^<l. 


de  d^e 
dz  af 


-K{e-e^y, 


6»<0. 


(9) 


Here  z^yt  and  ^  =  xl  Ij  .  The  cooling  conditions  are 
determined  by  the  constants 

K^lK^KFd,! p,E  and  51^  =(?; -7;)/(7; -T,) <0 •  (10) 


The  interface  resistance  enters  into  Eq.  (8)  as  a  combination 


R 


ll  E 


Kiz-E)- 


Rd,  ^ 


A 


(11) 


Here  7  is  the  current  transfer  length  [12]  that  determines  the 
distance  over  which  the  current  exchange  between  the 
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superconductor  and  stabilizer  takes  place.  Although  more 
complex,  Eqs.  (7-9)  share  common  features  with  the  standard 
Kardar-Parisi-Zhang  (KPZ)  equation  [14].  In  both  cases  the 
diffusion  competes  with  the  nonlinear  growth  termoc  {O'Y . 
The  KPZ  term  corresponds  to  the  heat  released  in  the  interface 
and  its  strength  is  proportional  to  the  interface  resistance. 

The  numerical  solutions  of  Eqs.  (7-9),  subject  to  periodic 
boundary  conditions,  were  obtained  for  an  initial  condition  in 
the  form  of  Gaussian  temperature  profile 

0{m  =  ^0  +(l  ^0 1  +^)exp{-f  /2^^} .  (12) 

All  results  presented  here  correspond  to  the  same  width  of  the 
initial  profile  5=V2.  The  value  of  the  dimensionless  ambient 
temperature  is  directly  related  to  the  transport  current  by  virtue 
of  the  relationship  (6): 

T(0) 

“  J 

where  jf^=J^{TQ)>J  is  the  critical  current  density  at  the 
operating  temperature  Tq. 


X 


Fig.  1 .  Sketch  of  the  conductor  and  two  examples  of  the  numerical  solutions 
of  Eqs.  (7-9).  Shown  are  the  evolving  temperature  profiles  0(x,t).  In  both 
cases  the  initial  condition  (in  the  foreground)  is  the  same,  Eq.  (12),  with 
00  =-l,  a  =1.1,  and  6  =1.4.  The  cooling  constant  k  =0.1  is  also  the  same.  In  (a) 
the  interface  resistance  is  negligible;  In  (b)  r  =  RIRQ  =1.5. 

III.  Resuets  and  Discussion 

A.  Speed  of  Normal  Zone  Propagation 
Figure  1  shows  two  examples,  in  the  form  of  2D  plots,  of 
the  solutions  0(x,  t)  of  Eqs.  (7-9),  subject  to  periodic  boundary 
conditions  and  the  initial  condition  given  by  Eq.  (12).  The 
normal  zone  evolves  from  the  initial  Gaussian  perturbation 
with  a  width  of  the  order  of  the  diffusion  length  (5  =1.4)  and 
the  peak  temperature  a  =  LI.  This  corresponds  to  the  peak  of 
the  real  temperature 


T^.=rHa-m-T,) 


(14) 


Fig.  2.  The  normal  zone  propagation  speed  as  a  function  of  interfacial 
resistance  obtained  from  the  solutions  of  Eqs.  (7-9)  similar  to  those  shown  in 
Fig.  1.  The  parameters  of  the  initial  condition  are  the  same  as  in  Fig.  1.  The 
adiabatic  case  corresponds  to  k=0.  The  value  of  0o=  -1  is  the  same  for  all 
cases. 


The  dimensionless  ambient  temperature  0o  =  -1,  that 
corresponds  to  the  transport  current  density  equal  to  50%  of 
the  maximum  critical  current:  j  =  .  Under  this  condition 

the  current  sharing  temperature  Ti,  defined  by  the  condition 
j^{T^)  =  lies  half-way  between  the  operating  temperature  Tq 

and  the  critical  temperature,  The 

maximum  temperature  of  the  normal  zone  is  determined  by 
the  cooling  constant, 

^NZ  .  (15) 

The  solutions  shown  in  Figs.  l(a,b)  correspond  toK=0.1, 

so  that: 


The  only  difference  between  the  solutions  in  Figs.  l(a,b)  is 
the  value  of  the  interfacial  resistance.  In  Fig.  1(a)  the 
resistance  r=0.  In  Fig.  1(b)  the  interfacial  resistance^  =  l.SR^ 


where  Rq  is  given  by  Eq.  (11).  The  propagation  speed 
determined  by  the  slope  dx/dt  at  the  level  of  0=1  (T=Tc) 
increases  drastically  with  increasing  interfacial  resistance. 

In  Fig.  2  the  propagation  speed  determined  from  the 
solutions  of  Eqs.  (7-9)  is  shown  as  a  function  of  the  interfacial 
resistance  for  three  levels  of  cooling.  The  NZP  speed  is  given 
by 


V  =  U,S{r,e„K), 
where  Ut  is  a  characteristic  speed, 

Kp,r 


(17) 


u,^i,r  =  (Drrr  = 


c\r-T,)d, 


(18) 

and  iS  is  a  dimensionless  speed  enhancement  factor  that 
depends  on  interface  resistance,  the  current  redundancy  0o, 
Eq.  (13),  and  cooling  conditions.  As  Fig.  2  shows,  the  NZP 
speed,  expressed  in  units  of  Ut,  increases  substantially  as  the 
resistance  of  the  interface  increases  above  the  threshold  Rq. 
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Fig.  3.  Stability  boundaries  as  functions  of  interfacial  resistance  for  three 
levels  of  cooling.  When  the  peak  temperature  of  the  initial  temperature  profile 
Tmax  given  by  the  parameter  a  in  Eqs.  (12,14)  exceeds  the  limit  shown  by  the 
curves,  the  normal  zone  starts  to  propagate. 

The  material  constants  typical  of  coated  conductors  are  as 
follows[13]:  Ky  ^  4-5  W/cm  K;  Q  «  1.7  jW  K;  pi  ^0.2x10'^ 
Q  cm,  and  t/i  =  40  pm  (copper  stabilizer).  K2^  7x10'^  W/cm 
K,  €2^  1.4  J/cm^  K,  and  ^2^  50  pm  (Hastelloy  substrate).  Let 
us  also  take  the  critical  temperature,  defined  by  the  condition 
Jc(Tc)  =0,  as  Tc  90  K  and  the  operating  temperature  To=65  K. 
As  an  example,  we  take  the  critical  current  density  Jc(To)  = 
200  A/cm,  and  the  transport  current  density  J=100  A/cm  (the 
corresponding  value  of  Oq  =  -1,  as  in  Fig.  2).  Then,  the  current 
sharing  temperature  Ti  ^  11  K.  These  parameters  determine 
the  increment  3.6  s'^  and  the  thermal  diffusivity  Dj  »1.4 
cm^/s,  which  gives  the  thermal  diffusion  length  ^  6  mm. 
Correspondingly,  the  characteristic  speed  Uj  =  Ijj  »  2.2  cm/s. 
Considering  that  for  low  interface  resistance  the  NZP  speed  is 
approximately  (  0.5-0.7)t/T,  see  Fig.  2,  these  estimates  agree 
well  with  the  results  of  the  experiments[4,5]  that  found  the 
NZP  speed  in  coated  conductors  of  the  order  of  1  cm/s. 

The  characteristic  interface  resistance  above  which  we  can 
expect  to  see  a  substantial  increase  in  propagation  speed  is 
Ro=pi(lTf/di^\'^^W^Q.  cm^.  The  typical  interfacial  resistance 
of  currently  manufactured  coated  conductors  is  of  the  order  of 
50  nQcm^[15].  Therefore,  the  interface  resistance  has  to  be 
increased  by  at  least  two  orders  of  magnitude  in  order  to 
observe  its  effect  on  NZP  speed  (see  Fig.  2). 

B.  Stability 

Stability  is  defined  as  an  ability  of  a  current-carrying 
conductor  to  dissipate  a  certain  amount  of  heat  and  return  to 
normal  operation  on  its  own.  The  NZP  speed  and  stability 
margins  usually  anti-correlate  in  superconducting  wires.  Our 
analysis  shows  that  the  effect  of  interfacial  resistance  is  no 
exception  from  this  rule.  With  increasing  resistance  the  NZP 
speed  increases,  but  the  ability  of  the  coated  conductor  to 
recover  from  a  perturbation  diminishes.  We  determined  the 
margins  of  stability  by  solving  Eqs.  (7-9)  with  the  initial 
condition  (12)  and  finding  a  peak  temperature,  determined  by 
the  constant  a  in  Eqs.(12,  14),  at  which  the  normal  zone  starts 


to  propagate.  The  width  of  the  perturbation  was  kept  constant, 
the  same  as  in  the  data  in  Figs.  1  and  2.  If  the  peak 
temperature  is  below  a  respective  curve  in  Fig.  3,  the  initial 
temperature  profile  dissipates  and  the  normal  zone  does  not 
propagate.  For  the  peak  temperatures  above  the  curve,  the 
normal  zone  expands. 

As  long  as  the  interfacial  resistance  R  is  below  the 
threshold  Rq  the  initial  temperature  perturbation  can  dissipate 
without  triggering  a  NZP  even  if  the  peak  temperature 
exceeds  the  current  sharing  temperature,  so  that  a  fraction  of 
the  transport  current  initially  flows  through  the  stabilizer. 
Even  when  the  peak  temperature  exceeds  the  critical 
temperature  the  conductor  remains  stable  if  the  cooling  rate 
determined  by  the  constant  k  is  sufficiently  large  (the  arrow  in 
Fig.  3).  With  increasing  resistance  the  ability  of  the  conductor 
to  recover  strongly  declines  and  for  R>  R^  a  perturbation 

with  a  peak  temperature  even  slightly  in  excess  of  T /  gives  rise 
to  a  NZP.  It  is  important  to  emphasize  that  at  i?  >  ^othe 

conductor  is  still  able  to  absorb  a  finite  amount  of  heat  and 
remain  stable,  but  only  as  long  as  no  current  is  diverted  into 
the  stabilizer  (within  our  model  it  means  T^ax^  Ti). 

IV.  Conclusions 

By  increasing  the  interfacial  resistance  by  about  two  orders 
of  magnitude  over  that  in  currently  manufactured  coated 
conductors  we  can  expect  to  see  a  substantial  increase  in  speed 
of  the  NZP.  A  side  effect  of  this  is  a  decrease  in  stability 
margins. 
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